Phonon-mediated electron pairing in graphene 
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The possibility of superconducting pairing of electrons in doped graphene due to in-plane and out- 
of-plane phonons is studied. Quadratic coupling of electrons with out-of-plane phonons is considered 
in details, taking into account both deformation potential and bond-stretch contributions. The order 
parameter of electron-electron pairing can have different structures due to four-component spinor 
character of electrons wave function. We consider s-wave pairing, diagonal on conduction and 
valence bands, but having arbitrary structure with respect to valley degree of freedom. The sign 
and magnitude of contribution of each phonon mode to effective electron-electron interaction turns 
out to depend on both the symmetry of phonon mode and the structure of the order parameter. 
Unconventional orbital-spin symmetry of the order parameter is found. 
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I. INTRODUCTION 

Low-energy dynamics of electrons in graphene, a two- 
dimensional form of carbon, is described by a two- 
dimensional Dirac-type equation for massless particles 
[l[. Such unusual electronic properties of graphene of- 
fer a possibility to study effectively ultrarelativistic elec- 
trons involved in condensed matter phenomena 13, |3|, 
and particularly in collective electron phenomena (j|. In 
the present paper, we consider qualitatively and estimate 
quantitatively Bardeen-Cooper-Schrieffer-like (BCS-like) 
[5| phonon-mediated pairing of electrons in graphene, 
taking into account the ultrarelativistic electron dynam- 
ics. The analogy exists between superconducting pair- 
ing of electrons in graphene and "color" -superconducting 
pairing in dense quark matter [6l|. 

As possible origins of electron pairing in graphene, 
phonon- and plasmon-mediated mechanisms [7| , electron 
correlations [a la] and highly anisotropic electron-electron 
scattering near van Hove singularity [lOj were proposed. 
Moreover, superconductivity can be induced in graphene 
due to proximity effect near superconducting contacts 
pdl . [l3 | . We study electron-electron pairing by in-plane 
optical phonons, represented by four modes with differ- 
ent symmetries, clearly seen in the Raman spectra of 
graphene (see, e.g., p^l - [l6| V Furthermore, we consider 
the quadratic interaction of graphene electrons with two 
modes of out-of-plane (flexural) phonons and study a 
possibility of electron pairing by these phonons (see also 
the qualitative study in [131 )• The properties of out-of- 
plane phonons and their interaction with electrons has a 
close relation to the formation of ripples in suspended 
graphene sheets [l^, [l3l and the the influence of rip- 
ples on electrons via effective gauge field [2^. The role 
of long- wavelength acoustic out-of-plane phonons in low- 
temperature transport of electrons in graphene was con- 
sidered in (211. 



Electron pairing in graphene, considered in the papers 
[22h27| within various models, can demonstrate various 
peculiarities, in particular, a multi-band character, when 
electrons from both conduction and valence bands are 
involved coherently into the pairing. Eliashberg multi- 
band equations for phonon-mediated pairing in graphene 
were derived and solved in [275 with neglecting details 
of electron-phonon interaction. Such details are taken 
into account in the present paper, where generally multi- 
band electron pairing is considered and the results for the 
superconducting gap in the one-band limit are presented. 

By means of matrix diagrammatic technique, we 
demonstrate that the effective interaction, induced by 
each phonon mode and entering the Eliashberg-type gap 
equations for electron-electron pairing, depends essen- 
tially on symmetry properties of this mode and on the 
structure of the electron Cooper pair condensate with re- 
spect to the valley degree of freedom (the analogue of chi- 
rality in graphene J28] ) . In result, the phonon mode can 
produce not only an effective attraction, but even effec- 
tive repulsion. Estimates of the coupling constants show 
that the quadratically-coupled out-of-plane phonons do 
not cause a pairing with any observable critical temper- 
atures, however the in-plane optical phonons can lead to 
the pairing in heavily doped graphene. 

The article is organized as follows. In Sec. 2 we for- 
mulate the Hamiltonian of graphene electrons, interact- 
ing with in-plane optical phonons. In Sec. 3, we derive 
from the first principles the Hamiltonian of electrons in- 
teraction with out-of-plane phonons. Electron-electron 
pairing due to in-plane and out-of-plane phonons is con- 
sidered in Sec. 4 and Sec. 5 respectively, and Sec. 6 is 
devoted to discussion and conclusions. 



II. HAMILTONIAN OF ELECTRONS AND 
IN-PLANE PHONONS 
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In this paragraph, we consider a Hamiltonian of inter- 
acting electrons and in-plane optical phonons in graphene 
as a base of subsequent consideration of the pairing. The 





FIG. 1: Graphene lattice as a combination of two triangular 
sublattices A and B with the period a — 2.46 A. The vectors 
dfe, {k — 1, 2, 3) and t; (/ = 1, . . . , 6) connect an atom of the 
sublattice A with its nearest and next-to-nearest neighbors 
respectively. 



FIG. 2: The reciprocal lattice of graphene in momentum 
space. The first and second Brillouin zones are encircled with 
thick solid and thin dashed lines respectively. The arrays of 
equivalent F, K and K' points are denoted by circles of dif- 
ferent colors. 



crystal lattice of graphene consists of two interpenetrat- 
ing triangle lattices A and B with a period a = 2.46 A 
(Fig. [Ij. We use the Hamiltonian of noninteracting 
graphene electrons in the tight-binding approximation 



Ho = -tJ2{atb,+b+a,} 



(1) 



where the sum (ij) is taken over pairs of nearest neigh- 
bors and t is the hopping amplitude; a^ and bj are 
destruction operators for electrons on the i-th and j- 
th sites of sublattices A and B. Performing in ^ 
the Fourier transform at — N^^^^^ e^^'^^ap, bj = 

7V-i/2^pe*P''^5p {N = 2S/a^V3 is the number of el- 
ementary cells in the crystal and S is the system area), 
we get 



Ho = -t^{a+bpgp + b+apg-p} , 



(2) 



where 
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and the vectors dk are shown in Fig. [TJ 

To describe low-energy electron dynamics in graphene, 
we focus on vicinities of the Dirac points K and 
K' = — K in momentum space (Fig. [5]). Introducing 
(analogously to [2y|) the four-component spinor opera- 
tor ^Jfp = (aK-i-p,&K+p,6K'+p,aK'+p)^, the"covariant" 
coordinates p° — {i/vp){d/dt), p^'"^ = Px,y, Pii = 
{p°, — p^, — p^} (here v-p — at^/?>/2 w 10^ m/s is the Fermi 
velocity), and the Dirac gamma- matrices in the Weyl rep- 
resentation 
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(3) 



where a — {(t^, ctj,, CTz}, we reduce ^ to the form: 



where ^p = ^p7'^. The operators ^p in the Heisenberg 
representation obey the Dirac-type equation: 



Pm7^*p = 0, m = 0,1,2. 



(5) 



The most general form of the Hamiltonian, describ- 
ing the linear coupling of graphene electrons to in-plane 
phonons, reads 



H, 



(lin) 
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(6) 
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where fj, enumerates the phonon modes, and gpq and F^ 
are corresponding coupling amplitudes and interaction 
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where 



is the phonon 



destruction operator. 

The in-plane phonon modes, most strongly coupled to 
electrons, are represented by Ai and Bi modes (we de- 
note them by fj. — I and 2 respectively) with the momen- 
tum q = ±K and energy wk ~ 0.150 eV, and by E2x 
and i?2y modes {fj, = 3 and 4) with q = F, wr ~ 0.196 eV 
I13l-ll6i. The interaction vertices for these modes 
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reflect their symmetry: scalar (/z = 1), pseudoscalar (/i = 
2) and pseudovector (/x = 3,4). The coupling constants 

ffpcj are weakly dependent on p, q: ^pq = g^, and can be 
related to the values (g|^ p)f, introduced in [IF 
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III. INTERACTION OF ELECTRONS WITH 
OUT-OF-PLANE PHONONS 

For out-of-plane phonon modes, we cannot restrict our- 
selves to phonon momenta close to F and ±K, as will be 
shown below. Therefore, we need a detailed microscopic 
model description of them. We start this description from 
the simplified version of the valence force Lagrangian 
[30l l31j , taking into account only the bond-bending term: 
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Here uf and u^ are out-of-plane displacements of carbon 
atoms from the sublattices A and B respectively, M and 
D are the carbon atom mass and the elasticity coefficient; 
Np{i) denotes the set of p-th order neighbors of the i-th 
atom. The Lagrange equations for atomic displacements, 
following from ([8]), are: 
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(9) 
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The displacements u^ ' can be decomposed over out- 
of-plane phonon modes with definite momenta q, covering 
the second Brillouin zone (see Fig. [2|), and branches a 
[a—1 for acoustical branch and 2 for optical one): 
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where e^^^ are polarizations of the phonon modes; $, 



qo- 



"--qfj 



-qo- 



is the phonon operator (in Heisenberg repre- 



sentation) and ujqa is its frequency. Substituting (|10| 
into (|9|), we obtain the system for determination of char- 
acteristics of phonon modes: 



M , ,2 u 
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where /iq = 12 + J2i=i ^**'*' (^^^ Fig. [T] for definition of 
tj). Solving (ITTt . we find: 
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where si = -1-1, S2 = —1. In principle, the summa- 
tion over q in (jlOp within the first Brillouin zone is suf- 
ficient, but it is convenient in our approach to use the 
second Brillouin zone, since the quantities (IT^ are de- 
fined unambiguously there. The long- wavelength asymp- 
totics of the frequencies ([T2|) are: uiqi « {a^ / 4) y^ D / M q^ , 
'^q2 ~ &^J D /M . Comparing cjq2 with the experimental 
data [13, we get ^D/M w 0.019 eV. 

Two mechanisms provide the main contribution to 
electron-phonon coupling in graphene: the deformation 
potential and the change of bond lengths [2l|, |3J|. In 
the simplest approximation, the deformation potential, 
acting on electron bound to specific carbon atom, origi- 
nates from the potentials of the nearest neighbors of this 
atom, approaching it (or moving away from it), and thus 
is determined by the sum of lengths of bonds, connecting 
this atom with its nearest neighbors. At the same time, 
the change of bond lengths modulate the hopping inte- 
grals in ([2]). In the limit of small out-of-plane displace- 
ments, the change of distance between i-th atom of the 
A sublattice and j'-th atom of the B sublattice is approx- 
imately {u^ — M^)^-\/3/2a. Therefore the Hamiltonian of 
quadratic electron-phonon coupling is: 



(»i> 



+C2{atbj+b+a,)], (13) 



where the multipliers of Ci and C2 correspond to the 
deformation potential and bond-stretch contributions; 
C2 = 7^3/20, 7 = dt/da « 6cV/A [H]. Taking the 
continuum long-wavelength limit of the deformation po- 
tential part of (IT^ and comparing it with [2l|, |33l, we 
find Ci = 51 /a2, where 51 sa 20 - 30 eV. 

Performing the Fourier transform for electron opera- 
tors in (fT3|) and using (fTO|) . we get 
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Here P and P' run over the first Brillouin zone, and 
q' = P' — P — q. For Li,L2 = A,B we have introduced 
the matrices: 
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Since electrons populate close vicinities of the momen- 
tums K and K' (Fig. [5]), we split the sums in ([Til) over 
P and P' among the valleys. Using the four-component 
spinor notations, we can rewrite (|14l) in the form: 

i^ir;f = ^E E *p'C'U'*p'J'<i'^'j'<i'-' (15) 

pp'q Qctct' 

where q' = Q + p' — p — q, the electron momentums p, p' 
are small and measured from the Dirac points, and the 
vector Q takes the values F and ±K. We can simplify 
cumbersome expressions for the vertices, assuming that 
they depend slowly on p and p'. Denoting Vq; ^ 



'00ci,aa' 



Vq^j, and using dH]), we get: 
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where 
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ELECTRON PAIRING BY IN-PLANE 
PHONONS 
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FIG. 3: Diagrammatic representation of self-consistent gap 
equations: (a) for linear electron-phonon coupling (|19p . (b) 
for quadratic electron-phonon coupling (|3ip . 



Employing the standard diagrammatic technique for 
the system of graphene electrons with the Hamiltonian 
(HI, interacting with in-plane phonons via ([6]), we get 
the following set of Gor'kov-type equations for the Green 
functions (J17p . describing the pairing in the mean- field 
approximation: 



G„- ip) = gI? (p) - ^ E 9lD,{p - P')G^^ ip) 

xF|;)G,,.(p')r('=^Gfc,(p), (18) 



where F^'^ = r^, F|f^ = 
conjugated vertex and D^{q) 

phonon Green function. Introducing the anomalous self- 
energies 



G^^FjG is the charge- 
) is the 
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for i ^ j, we rewrite ()18p in the form 



(i) 



(19) 



To describe pairing, we introduce the following set of 
matrix (4 x 4) Green functions in Matsubara representa- 
tion (similarly to [34 



G„(p,t) 



,(1) 



-{T¥;Hr)^p\0)), 



(2), 



(17) 



C^_ 



where *^'^(t) = *p(r), *^-^^(t) = *cp(t) 
is the charge-conjugated spinor and C — iy^^'^ is the 
charge-conjugation matrix. The diagonal elements of 
([T7[) . Gil and G22, are the Green functions of parti- 
cles and "antiparticles" , while G12 and G21 describe a 
Cooper pair condensate. The Green functions of free 
particles are: G^°^(p) = [7°(po + M)-'yF7P]"\ G22^(p) = 
[7"(Po-m)-^'f7P]"\ where p = {po = z7rT(2n-h l),p} 
and /I is the chemical potential in graphene, measured 
from the Dirac points. 



^-1 _^(o)-i . r^(o)A 
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1^12 — (^11 ^12'J22, (J21 — <-^22 ^21'Jll- 



A2lGj^j^ A12, 



(20) 



The diagrammatic representation of (|19p is shown in 

Fig.ETa). rn i-in 

Similarly to [23 . |3J-[3a|, we employ the following 
method to solve the matrix equations ([T ^ ~ (|20p : 1) we 
assume a certain form of the order parameter A^^- , which 
has a definite matrix structure and is parameterized by 
a set of variables; 2) with this A^, we solve the first 
pair of the Gor'kov equations ([^ and find the inverted 
normal Green functions G~^\ 3) inverting G~^ , we sub- 
stitute Gii into the second pair of the equations (|20p . 
finding the anomalous Green functions Gij ; 4) substitut- 
ing Gij into ([ro[) , we get a closed system of equations for 



Aab) 



variables, parameterizing Aij. The matrix structures of where C±(p) 
Gij, found from (PO]) . should correspond to that of the 
initially assumed Aij via (fT9|) . and generally this occurs 
only for the certain forms of the order parameter A,j. 
Possible matrix structures of the order parameter, de- 
scribing the pairing of relativistic elementary particles or 
electrons in graphene in various models, were discussed 
in [22,, 25, 26, 34 38,]. 

For phonon-mediated electron pairing, we assume the 
simplest form of the order parameter: we suppose that 
the pairing is s-wave and diagonal with respect to con- 
duction and valence bands. In the relativistic-like ap- 
proach to the electron dynamics (O, the states of elec- 
tron in the valleys ±K are the states with the "chi- 
rality" quantum numbers ±1 respectively. Moreover, 
the electron states in conduction band (positive-energy 
"particles") have the equal chirality and helicity, while 
the electron states in valence band (negative-energy "an- 
tiparticles" ) have the opposite chirality and helicity. The 
projection operators on the states with definite chirality 
(c = ±1) and helicity {h = ±1) are [H: 



— itrlpl — M Emd E±{p) = 

^± (p) + A j_ {p) are the excitation energies for bare 

electrons and Bogolyubov quasiparticles in conduction 
and valence bands. Using ([^ and ([^ in P^ . then 
multiplying the both parts of (ITOl) by 7'a(p)e^''^'^ from 
the left and taking a trace, we derive the system of 
self-consistent equations for two gaps A±(p): 






^l^g^^D,{p-p')-^ 



A^pQ 



pV /3=± 



Elip') 



xii;^^(p,p';v), (24) 
where the following angular factors are introduced: 
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(25) 
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*7 7 7 7 J P 



p/|p|. The following oper- 



ators make projections on the conduction (a 
valence (a ~ —1) bands: 



-fl) and 



The contact character of electron interaction with the 
optical phonons allows to perform angle integration 
over p and p' in (P5)) . which leads to independence of 
_R^„(p,p';v) = Rfj.i'v) over a and /3. Calculating the 
angular factors ((25|) using ^, ((2T|) and summing over 
degenerate modes, Rk = Ri + R2 and Rr = R3 + R4, 
for K- and F-phonons respectively, we get: 



^q(p) = Vc^a-hip) + V-c-a-hip) 



1 -I- a7°7p 



(21) 



i?K = — cos^ V + -vl sin^ V, 
i?P = - cos^ V + (-Vj + v| + V3) sin^ v. 



We assume an arbitrary structure of the order param- 
eter with respect to valley degree of freedom, represented 
by (2 X 2) unitary matrix in the S£ace of the valleys 
{K, K'} (the same was supposed in [23j, where pairing of 
Zeeman-split electrons and holes in graphene was stud- 
ied). The operators Ti = 7^, Tz = 7^75 and T3 = i-/^ 
allows us to construct the expression for Aij explicitly. 
These operators perform rotations of the order parame- 
ter in the space of valley degree of freedom and, as noted 
in [201, obey the algebra of the Pauli matrices. Thus, the 
general form of the order parameter A21, corresponding 
to the band-diagonal pairing with the arbitrary valley 
structure, is 



(26) 



The system of gap equations (j24p can be rewritten in 
the form of two-band gap equations [25l-[27| 



Aa(p)-^EE^^ff(P-P» 
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p' /3=± ''" '^P^ 

but with the effective interaction, dependent on v: 
T4ff((?;v)=^52^^(g)i?^(v). 



A2l(p) 



^JvT 



A+(p)n(p)+A_(p)P_(p) , (22) 



where A+(ip) are the gaps in conduction and valence 
bands (see [2^, [2g|), and the three-dimensional vector 
V parameterizes the valley structure of the order param- 
eter. 

Substituting (|22|) into (pO)) and taking into account, 
that A12 = 7"A+ 7°, we find: 

r^ ( \ S^ Po+Ca(p) -5 ^M 

V^ Aa(p) 

<^2l{p) = > -2 ^^TTT 



7°e-T^^(p)^o^ (23) 



(27) 



(28) 



The effective interaction ([28l) is a linear combination of 
interactions due to separate phonon modes with the co- 
efficients i?^, which can take values in the range — 1 < 
R^ < 1 depending on v. Therefore, the in-plane optical 
phonons in graphene can induce not only attraction be- 
tween electrons {R^ > 0), but even repulsion (i?^ < 0). 

Consider now the limiting cases of the valley structure 
of the order parameter. At v = 0, the valley part of the 
order parameter is the unit (2 x 2) matrix in the valley 
space; in this case, all -R^ = —1, and all the phonon 
modes lead to effective repulsion. When v — {tt/2, 0, 0}, 
the valley part of the order parameter is ct^ , and i?K = 1 
(effective attraction), Rr = — 1 (repulsion). Finally, at 
V = {0, 7r/2,0} and v = {0, 0,7r/2}, the valley part is 
ax and ay respectively (valley-off diagonal pairing) , and 
^K = (mutual cancelation of contributions from scalar 
/i = 1 and pseudoscalar ^ — 2 phonon modes), i?r = 1 
(attraction). 



The maximal value of effective electron-phonon cou- 
pling constant can be reached in highly doped graphene 
due to large density of states at the Fermi level N = 
/i/27rwp. At /i > oj^ the pairing can be treated as one- 
band (i.e. involving only the conduction band). In [27|, 
the equations of the type of ([27]) were considered and cor- 
responding Eliashberg equations were derived and solved 
both in the one-band regime and in the limit of small 
graphene doping. The Eliashberg function a^F, corre- 
sponding to (pS)) . is 



coupling results in the loop of two phonon lines con- 
necting two interacting electrons [l3, [2l[ , instead of one 
phonon line for linear coupling, as shown in Fig. EJb). 
The analogue of self-consistent gap equations ([T^ for the 
electron-phonon Hamiltonian (llSp reads: 
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(^-iy(Q)T ^^ jg ^j^g charge- 



and the equation for the gap in the conduction band A — 
A+(pf = m/'^f,'^ = 0) at T = reads (see [22 
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where wg 



duj 



Vu;2 _ 1^2 



dv 



al{v)F{v) 



(30) 



is a cutoff frequency for the gap 



A+(pf,w). Solving (pO]) with taking into account ([29|) 
in the limit A <C ujq^uj^, we obtain the estimate of the 
gap at T = 0: 



where V^.^^„„, - ^ ^-p',-pq,aa' 
conjugated vertex. The charge conjugation does not 
change the contributions to the vertices p^ from bond- 
stretching, but changes the sign of the deformation po- 
tential contributions. Note, that the summation over q 
in pil) is performed over the whole second Brillouin zone. 
Performing a summation over go = 27rOT in the phonon 
loop in (1311) , we get an analogue of phonon Green func- 
tion, but with the sum of the frequencies of two phonons: 
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Further calculations are similar to that in derivation of 
the Eliashberg equations (|29p - ([5n)) . We consider one- 
with the partial A^ = 2NglR^{w)/uj^ and total A = band pairing again, occurring at // > cjq^ - 0.1 eV. 



'^ Xfj, coupling constants introduced 

The valley structure of the order parameter v will ad- 
just itself to achieve the ground state with the lowest 
energy, corresponding to the maximal possible value of 
A. When (^k/'^k) > 2(gp/(u;r), the preferable pairing 
structure is valley-diagonal: v = {7r/2,0, 0}; in the op- 
posite case, when (g'^/uj-K) < 2((7p/ajr), the pairing is 
valley-off diagonal: v — {0, {n/2) cos(p, {n/2) simp}. 

Taking the values of the coupling constants from [13| , 

we get: g^/uJK w 12.02 eV • A^ g^ujr ~ 3.75 eV • A^ 
so the valley-diagonal pairing with v = {7r/2, 0, 0}, when 
F-phonons compete with the K-phonons, is preferable 
(note, that this relation can revert at large dielectric con- 
stant of surrounding medium, since g^ is highly renor- 
malized to higher values due to Coulomb interaction 
[l5|). The coupling constant A ~ 0.065 x ^ [eV] can 
provide any noticeable pairing only at heavy chemical 
doping of graphene with /x w 1.5 — 2eV. The earlier 
estimates [Ij, [23| of the coupling constants for electron 
pairing in graphene, mediated by optical phonons, give 
similar results. 



V. ELECTRON PAIRING BY OUT-OF-PLANE 
PHONONS 

The consideration of electron pairing by out-of-plane 
phonons is based on the same approach, as in the previ- 
ous paragraph, but with using the interaction Hamilto- 
nian (J15p instead of ©. The quadratic electron-phonon 



The Eliashberg function, corresponding to out-of-plane 
phonons, and entering the equation analogous to (j30p . is 



M 
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-V 
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^UgW.^Uy-^ 



qcrcr 



(32) 



Here we neglected the momentums p and p' in the ver- 
tices and phonon frequencies, as in (1161) . and performed 
angle integration. 

It is convenient to rewrite (1301) in the form: 



1 = 



duj 



v^ 



Z{lj]v), 



Z{uj; v) = du 






(33) 



Summation over q within the second Brillouin zone in 
(15^ with taking into account ([T5| results in some ef- 
fective matrix operation with the expression 7*'e*^'^7'^, 
enclosed between the vertices. A number of terms un- 
der the trace in p2l) vanish upon summation, and the 
remaining nonzero terms allow to write p3p as a sum 
of contributions, corresponding to deformation potential 
(/i = 1,2) and bond-stretching (/i = 3,4) and, on the 
other side, to the two-phonon processes, leaving the elec- 
tron in its initial valley (/x = 1, 3) and flipping it into the 




tions Z^J^ of x 



Lo/^D/M, plotted in Fig. gl The func- 
tions Zi and Z3 have logarithmic singularities at oj = 
due to contributions of acoustical branches, indicating 
long-range character of electron-electron interaction by 
out-of-plane phonons (which was noted in [l7[), but gen- 
erally all Z^ are of the order of unity. Therefore, the 
summary function Z[uj), playing the role of the coupling 
constant in (p3| is of the order of 10"'^ at v = {0, 0, 7r/2}, 
and cannot provide any observable pairing. 

Note that consideration of electron pairin g in . graphene 
due to out-of-plane phonons, performed in 17|, neglects 
details of electron-phonon interaction and includes only 
acoustical phonon branch, but provides the order of mag- 
nitude for the coupling constant, close to that in our 
work. 



FIG. 4: The dimensionless functions Zfj,{x), corresponding to 
Zfj,{Lj) (|34p . Solid line: fi = 1, dashed line: fi = 2, dotted line: 
/i = 3, dash-dotted line: fi — 4. 



opposite valley (/i = 2,4): Z(w; v) = J2l=i ^m(^)^m(v)- 
Explicitly, we have 

i?i = 1, i?2 = (-V2 +V3)sin^i;, 

n 2 / -2 -2 -2\-2 

R3 = — COS V — (V]^ ^ ^2 — V3) sm V, 
i?4 = — cos^ V + -vl sin^ V, 
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(34) 



The structure of the order parameter, most fa- 
vorable with respect to all branches ^ = 1 , . . . , 4 
jointly, is V = {0,0, 7r/2} (valley- antidiagonal and 
antisymmetric pairing). In this case, i?i,2,3 = 1, 
R4 = 0. The functions ([M|) can be reduced, 
Zi,2(^) = {3UgyWaHl^/^D^/^)Z,^2{x), Zs^u:) = 
{^M^"^ li2M^I'^D'^l'^)Z-i,^{x), to the dimensionless func- 



VI. DISCUSSION 

In addition to the previous sections, we put few gen- 
eral remarks. First, the issue of symmetry properties 
of the order parameter with taking into account elec- 
tron spins is worth of discussion. The spin projection 
indices s, s' can be assigned to the anomalous Green 
functions, Gi'}^(p,T) = -(T*cp«(t)'^ps'(0)), to give 
the following condition of antisymmetry: G^^, (p,0) = 

C~^[G^,g'(— p, 0)]-^C. Applying it to the order param- 
eter of the form ([22|) . we get that when the pairing is 
valley-diagonal (i.e. v ex {1,0,0}) or valley- antidiagonal 
and valley-symmetric (v ex {0,1,0}), the order param- 
eter must have combined spatial-spin symmetry (spin- 
triplet s-wave pairing) , unlike the usual electron-electron 
pairing [5|, which is jointly antisymmetric in the space 
and spin (spin-singlet s- or d-wave pairing, or spin- 
triplet p-wave pairing). Such a peculiarity, as can be 
shown, is a consequence of additional "hidden" antisym- 
metry of the order parameter by sublattices (the similar 
unconventional symmetry of two-electron wave function 
in graphene was noted in [33] )■ Conversely, when the 
pairing is valley-antidiagonal and valley-antisymmetric 
(v ex {0,0,1}), the order parameter must be antisym- 
metric jointly in space and by spins, as in conventional 
superconductors . 

Another important point is a role of Coulomb repulsion 
of electrons, which can be added to the self-consistency 
equations (HH) with the vertices F^^) = -F^^) = 7". The 
corresponding angular factor (j26p does not depend on v 
and always provide the effective repulsion in the frame- 
work of the band-diagonal pairing ([2^ . However, cer- 
tain forms of the order parameter are possible, support- 
ing the electron-electron pairing by Coulomb interaction 
in graphene: for example, the "vector" order parameter 
in [22[ or the resonating valence bond order parameter 
[si |9| ; all of them can be described using the formalism 
of matrix Green functions (ITTt . 

When graphene is heavily doped (/i > 0.5 eV) by im- 
purities (see, e.g., [40|), several additional factors should 
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be taken into account when considering electron pair- 
ing. These are, in particular, an influence of impurities 
on the condensate [41| , formation of energy bands of the 
deposed atoms (similarly to that in graphite intercala- 
tion compounds 1421 ) and possible structural reconstruc- 
tion of graphene [40|. Moreover, the trigonal warping of 
the Fermi surface in graphene at high doping [l| should 
promote the valley- antidiagonal pairing of electrons with 
opposite momenta. 

In conclusion, we have considered phonon-mediated 
electron-electron pairing in graphene taking into account 
both the details of electron-phonon interaction resolved 
by sublattices and valleys, and a possibility of different 
structures of the order parameter. We assumed an s- 
wave pairing, diagonal with respect to electron bands, 
and shown, that conditions of this pairing depend on a 
structure of the order parameter with respect to the val- 
ley degree of freedom. Contribution of phonon modes 
in graphene to the effective electron-electron interaction 
can be attractive, repulsive, or can even vanish, depend- 
ing on both the mode symmetry and the valley structure 
of the order parameter. The orbital-spin part of the order 
parameter can be symmetric in some cases. 

We have also considered the quadratic coupling of 
electrons to out-of-plane phonon modes. Hamiltonian 
of electron-phonon interaction is derived taking into 
account the contributions from deformation potential 



and from change of bond lengths. Quadratic character 
of electron-phonon coupling results in unusual phonon- 
mediated electron-electron interaction, represented by 
the loop consisting of two phonon lines. Integration on 
the inner momentum in this loop should be performed 
over the whole Brillouin zone, in contrast to the situa- 
tion with the linear electron-phonon coupling, where the 
phonon momentums are either very small, or connect two 
electron valleys. The effective action of the out-of-plane 
phonons on electrons is also dependent on the structure 
of the order parameter. 

The study of phonon-mediated electron pairing in 
graphene, presented in this paper, not only allows to ex- 
tend the analogy between electrons in graphene and rela- 
tivistic elementary particles, introducing new kind of in- 
teractions in the model (scalar, pseudoscalar, pseudovec- 
tor phonons etc.), but would provide better understand- 
ing of BCS-like pairing phenomena in unconventional sys- 
tems. 
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